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Abstract 

The quark-loop contribution to the p° — u> mixing self-energy function is cal- 
culated using a phenomenologically successful QCD-based model field theory 
in which the p° and to mesons are composite qq bound states. In this calcula- 
tion the dressed quark propagator, obtained from a model Dyson-Schwinger 
equation, is confining. In contrast to previous studies, the meson-qq ver- 
tex functions are characterised by a strength and range determined by the 
dynamics of the model; and the calculated off-mass-shell behaviour of the 
mixing amplitude includes the contribution from the calculated diagonal me- 
son self-energies. The mixing amplitude is shown to be very sensitive to the 
small isovector component of dynamical chiral symmetry breaking. The space- 
like quark-loop mixing-amplitude generates an insignificant charge symmetry 
breaking nuclear force. 

PACS numbers: 12.38.Lg, 13.75.Cs, 11. 10. St, 12.38.Aw 
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1. Introduction. The on-mass-shell mixing of p° and u mesons can be represented via 
a 2 x 2 matrix, Ai 2 , with diagonal elements m 2 and and both off-diagonal elements equal 
to the mixing amplitude IT. The established value of IT, deduced from e + e~ — > 2tt data, 
is —4520 ± 600 MeV 2 JTJ. Recently [@ the systematics of electromagnetic and semi-strong 
mass splittings of the hadrons have been used to deduce a universal strength of —5000 MeV 2 
for isovector mixing; i.e., for both p° — u and 7r — rj. The contributions that this mixing 
makes to the various charge symmetry breaking (CSB) phenomena in nuclear physics have 
been much studied [|l],[3|]. The central element is the corresponding CSB contribution to the 
nucleon-nucleon potential, which is obtained in the boson exchange approach by using 



(k 2 — m 2 ) (k 2 — m 2 ) 

as the off-diagonal propagator. Without evidence of significant momentum dependence, II 
has been taken as a constant in nuclear CSB applications and the resulting phenomenology 
has been very successful. In a number of circumstances, the p°—u CSB nuclear potential 
thus obtained has been seen to dominate over other mechanisms, such as electromagnetic 
effects, the mass difference of neutrons and protons and 7r —r] mixing. A particularly striking 
example is the IUCF measurement [§[] of the analyzing power difference in np scattering at 
180 MeV. 

In QCD, the mixing arises because of the small u-d current-quark-mass difference, which 
induces a difference between the propagation of the uu and dd components, and one expects 
a momentum dependent off-diagonal self-energy. Knowledge of the off-mass-shell behaviour 
is needed to extrapolate to spacelike momenta appropriate to meson exchange between nu- 
cleons. The simplest possible contribution to U(k 2 ) is a single quark loop and the initial 
investigation employed free constituent quark propagators and phenomenological meson-g- 
q- vertex form factors to obtain a closed form expression ||. The relevance of the result is 
clouded somewhat by the presence of a spurious quark production threshold, inconsistent 
with the notion of quark confinement. To avoid this unphysical feature, which arises near the 
meson mass-shell, a subsequent work || employed a simple, phenomenological form of dy- 
namical quark propagator that implements confinement through the absence of a mass-shell 
pole. Nevertheless the central result remained the same: for spacelike momenta appropriate 
to vector exchange between nucleons, the mixing is so much weaker than the mass-shell 
value that the corresponding CSB nuclear potential is negligible. A recent QCD sum rule 
investigation |§ found a somewhat stronger momentum dependence, with U(k 2 ) reaching 
the negative of the mass-shell value at fc 2 ~ 0. In contrast, the momentum dependence 
of the 7T° — rj mixing amplitude calculated from chiral perturbation theory appears to be 
considerably weaker J/J. 

There are a number of issues that suggest the quark loop calculations may not yet be 
reliable for this process. An important consideration for p° — uj mixing is that it is driven 
by the small isovector component of the quark propagator, which in turn is driven by the 
isovector current mass m u — rrid. A direct parametrisation of a propagator model is usually 
constrained by a fit to quantities related to dynamical chiral symmetry breaking (DCSB), 
such as fn, < qq > and m n . These are isoscalar constraints and are not sufficient here. The 
explicit chiral symmetry breaking (ECSB) generated by current masses is weaker but its role 
is magnified in the isovector propagator and some dynamical guidance is to be preferred. 
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The low momentum behaviour of the dressed quark propagator can be calculated using 
a dynamical model Dyson-Schwinger equation (DSE). However, typical DSE studies are 
carried out numerically in Euclidean space; i.e., at spacelike momenta, where maximal use 
can be made of the constraints from perturbation theory and the renormalisation group in 



QCD P,l0fl. Such numerical studies do not provide readily accessible information about the 
behaviour at real, time-like quark momenta nor the domain of complex momenta sampled 
in a Euclidean treatment of meson propagators at time-like meson momenta. 

The above issues are addressed in this work through use of a recently developed solution 
|TT| of a model DSE that yields the quark propagator in closed form as an entire function in 
the complex momentum plane. This structure can be interpreted as representing confined 
quarks [§,|T3 an d provides an unambiguous representation over the complex momentum 
domain required by the quark loop. The dependence upon the current mass is explicit in 
the model and is produced by the DSE dynamics. 

We calculate the mixed p° — u self-energy function produced by the vector meson sector 
of a QCD-based model field theory pHH, the Global Colour-symmetry Model (GCM) [for 
a review see ||15|| 1, which formalises the coupled Dyson-Schwinger-Bethe-Salpeter equation 
approach to QCD phenomenology and allows several internally consistent features to 
be implemented for the first time. The off-mass-shell behaviour we obtain for the mixing 
amplitude includes the contribution from the off-mass-shell structure of the composite p and 
uj propagators. In addition, the strength and range parameters used to model the meson-gg 
Bethe-Salpeter amplitudes are determined by the meson mass-shell dynamics within the 
GCM. This eliminates the need for several assumptions made in previous studies. A short 
account of some of this work has recently been presented [JT6 . 



2. p°-uj Mixing Amplitude. The description of the p° and u mesons as q-q composites 
is obtained from a phenomenologically successful, QCD-based model field theory |I3,I3J, 
which is defined by the action: 

S[q, q] = / d 4 x q{x) (7 ■ d x + M) q(x) + \ [ d 4 xd A y j a (x)g 2 D(x - y)f(y), (2) 



where M is the current-quark-mass matrix, j®(x) = q{x)^^^q{x) is the quark current, and 
D(x — y) is a model effective-gluon-propagator in a Feynman-like gauge. Here we consider 
only u and d flavors and use a Euclidean metric throughout such that a ■ b = a^b^ and 

{7m> lv) = 2 V> with 7m = 7j- 

After bosonisation and expansion about the classical vacuum [by which is meant the 

minimum of the action], the tree-level effective action for the mesons is |f["7H 

00 / j\n 

S = TrJ2 [^0(^75^ • 7? + hvUv + ijvT • Pv H )]" 



„=2 n 



2g2D(x _ y ] — , (3) 

where the meson fields are bilocal fields corresponding to bilocal combinations of quark 
fields; for example, oo^(x,y) ~ q{y)i^ '^q{x) . The dressed quark propagator, G , that ap- 
pears here has a self-energy determined by the classical vacuum configurations of the bilo- 
cal fields. We use a momentum representation in which q is conjugate to x — y and 
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P is conjugate to (x + y)/2. The vector bilocal fields are represented by, for exam- 
ple, l^^q-P) = T(q;P)^T^(P)co u (P), where T^P) = 5^ - P^PJ P 2 and T(q;P), a 
Lorentz-scalar function, is the dominant Bethe-Salpeter amplitude in the vector qq channel 
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We are interested in tree-level coupling in the p° — uj sector. In a matrix notation where 
Vfj, denotes the column of transverse fields the GCM action up to second order in 

the fields is 

s[pM = \ J *?Lv?{-p)[a£(p) + n M ,(P)]K(P), (4) 

where A -1 is diagonal and IT is the off-diagonal p° — uj self-energy. The diagonal inverse 
propagator is 

A^(P) = / -^tr [G ( q -)ij»fGa(q + )i lv f) T 2 (q; P) + 95, u J ^fejy. (5) 

with tr denoting a trace over spin, flavor and color. Here / are the relevant flavor factors 
[1 for the uj, and f for the p\ and q± = q ± P/2, where P is the meson momentum and q is 
the loop momentum associated with the qq substructure. The mixed self-energy, II M j,(P), is 
given by an expression identical to the first term of (||) except that one of the flavor factors 
is T3 and the other is 1, thus it is nonzero only if Go has an isovector component. 

At tree level the diagonal terms are the same for p° and uo and the general form of the 
transverse component is 

A-J(P 2 ) = T^(P) (P 2 + My) Z(P 2 ), (6) 

where the position of the zero defines the physical mass My. The meson fields may be 
rescaled so that at least the on-mass-shell value of \fZ is absorbed into the fields to pro- 



duce a physical normalization. We elect to absorb yZ(P 2 ) into the fields so that the 
diagonal propagators have the standard form characterised only by the physical mass. The 
corresponding mixed self-energy U^ U (P)/Z(P 2 ) then contains the off-mass- shell structure 
appropriate for use in meson exchange applications that use the standard point-meson form 
for diagonal propagators. The net result of this is that the renormalised version of (|5|) is 
obtained via the replacement: 



The resulting vertex function, F(q;P), contains the strength of the meson-gg coupling as 
determined by the model. On the mass-shell this is equivalent to imposing the standard 
physical normalisation for a Bethe-Salpeter amplitude (TP] . 
The corresponding mixed self-energy 

IV(P) = / -0^tr [Ga(q-)i lfi T 3 G (q + )i lv } F 2 (q; P) (8) 

has a transverse component, 11^ (P) = T^(P)Ut(P 2 ), which yields the mixing amplitude 
as the following difference of quark loop integrals: 
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IT T (P 2 ) = W T {P 2 ) - IT|(P 2 ). (9) 

With the quark propagator in the form Go(p) = —27 • p<Jv(p 2 ) + &s(p 2 )i the required quark 
loop integrals are 

W T (P 2 ) = (10) 
12 



(27T) 



/ dssF 2 (s; P 2 ) jf_ i - u 2 )^[(l (1 + 2u 2 ) - |P 2 )<(+)<(-) + <tJ(+KR 



where s = q 2 , u = q ■ P and cr(±) denotes <j(q±), where q± = s + P 2 /4 ± uy/sP 2 . For time- 
like meson momentum, where P 2 < 0, yfp 2 — *> i\/—P 2 and cr(— ) = cr*(+). By assumption 
the vertex function, P, is independent of the angle u. 

3. Quark propagator. Because the momentum region of interest [P 2 > — m 2 ] includes 
a time-like sector, the loop integral fllPf ) requires knowledge of the propagator amplitudes 
<7 s (z) and <7 W (^) in the domain of the complex z = x + iy plane enclosed by the parabola 
\y\ = vn>u \Jx + in 2 1 4 with x > — m 2 /4. A fit to soft chiral quantities [such as (qq), f^, m^, 
r n ] does not probe a propagator model over this domain. [The corresponding quark loop 
appropriate for pion physics also requires knowledge of the quark propagator in the complex 
plane but in a region whose area is reduced by a factor of {m^/m^) 2 ~ 25 and hence the 
chiral limit (m T — > 0) captures the dominant physics.] 

The quark propagator given by Ref. []TTJ| provides a behaviour in the complex domain 



that is accountable to DSE dynamics. It is conveniently expressed in terms of dimensionless 
amplitudes a s (x) = \a s (q 2 ) and a v (x) = X 2 a v (q 2 ), where x = q 2 /\ 2 and A is the momentum 
scale. The amplitudes a are also explicit functions of the dimensionless current-quark mass 
fh = m/X. This DSE model employs a momentum-space delta function for the gluon 
propagator, with a strength set by A 2 , and a dressed quark-gluon vertex p0| , which is 



compatible with the Ward-Takahashi identity. The amplitudes a are entire functions in 
the complex momentum plane with an essential singularity at timelike infinity. The exact 
expressions are a little complicated when the current mass is included. We have verified 
that for typical current masses [~ 10 MeV], the following simple forms 



TO 

a s (x) = c(m)e- 2x + -(1 - e~ 2x ), (11) 



x 



ovtx) = 6 ~ 2X ^ 2X) - m c(m)e- 2 *, (12) 

are accurate representations [to order 10~ 4 ] on the domain of complex momenta encountered 
in the quark loop. 

These forms also accurately represent the leading behaviour of the quark propagator 
at large spacelike q 2 in the model and in QCD. The constant c(m), which represents the 
strength of DCSB, is not determined by the model DSE solution; the reason being that 
the model of Ref. |l ij has only two scales [momentum A and mass m]. A third scale is 
necessary to fix c(rh). In a more elaborate model, such as that of Ref. |TIJ, this is provided 



by Aqcd, which is introduced by the inclusion of the known large spacelike-g 2 behaviour of 
the gluon propagator. However, a numerical solution is then necessary and the extension 
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into the complex momentum plane is not unique. We consider ([11]) and (|12|), with c(m) a 
parameter to be determined below, to provide a realistic model propagator for a confined 
quark in QCD over the region of the complex plane sampled in ([TDD • 



To fix the values of the parameters A and Co = [c u + q]/2 we follow Ref. |H|] and 
use (]TTJ) and flTjP to fit (gg), J^, 7V and the 7r-7r scattering lengths: a[j, a^, a{, a°. With 
(m u + md)/2 = 7.5 MeV this leads to best fit values of A = 0.893 GeV and Co = 0.523, which 
yield (uu) = (dd) = (183 MeV) 3 , f n = 90 MeV, r T = 0.56 fm and a mean deviation of 19% 
from the experimental scattering lengths. 

4. Bethe-Salpeter amplitude. The dominant vector meson Bethe-Salpeter amplitude, 
F(q; P), also appears in ([T~0[) . This may be calculated in the GCM using a variational form of 
the off-mass-shell generalised-ladder-approximation Bethe-Salpeter eigenvalue equation WE 



and a simple, reasonable approximation for the shape on the mass-shell is provided by the P- 
independent form T(q\ P) = exp(— q 2 /a 2 ), with a = 0.62 GeV. To simplify the computations 
here the renormalisation function, Z(P 2 ), was obtained from ([5]) after the second [constant] 
term is set to produce the vector mass 783 MeV. The normalised amplitude F(q; P) is then 
obtained from (0). 

5. Isovector component of the quark propagator and p°-uj mixing. It only re- 
mains now to set the isovector component of the quark propagator, which is a measure of the 
isovector character of DCSB. With the representation Gq 1 (q 2 ) = ■ qA(m, q 2 ) + B(m, q 2 ), 
one commonly used assumption is that A(m, q 2 ) ~ A(0,q 2 ) and B(m,q 2 ) « B(0, q 2 ) + m. 
This is the procedure followed in the one other calculation of p° — uo mixing that used a 
dynamical quark propagator ||. It is convenient because only the chiral limit amplitudes 
need to be modelled. With the m = forms of A and B extracted from flTlD and (|T2|), this 



yields the prediction for the p° — uj mixing amplitude shown as the dotted line in Fig. 1. 
[Throughout this work we use — m u = 5 MeV.] 

The momentum dependence obtained is quite similar to that found in earlier work . 
The experimental mass-shell mixing strength was produced in the work of Ref. 0] with 
use of phenomenological meson-quark coupling strengths g p = 4 and = 5 deduced by 
simple scaling arguments from empirical meson-nucleon coupling constants. In the present 
approach, however, the vertex strength is an outcome of the model, and the comparable 
measure from (|7|) is g p = g w = F(q = 0; P 2 = —My) which produces the larger value 16.4. 
Nevertheless, we obtain only 50% of the experimental mixing value. The quark vertex 
alone is not an accurate guide to the net coupling strength between hadrons because it 
does not account for the particular way the hadron dynamics samples the quark propagator 
amplitudes, especially the wave function renormalisation function. We feel the mass-shell 
condition used here provides the internally consistent vertex normalisation for the meson 
model. 

The above estimate for the isovector propagator can be improved through the use of 
solutions from a realistic model DSE to directly determine the isovector component of c(m) 
in (|TTJ) and (|T2"D . [We note that in realistic DSE studies A{m, 0) decreases with increasing 



m and B(m, 0) is badly underestimated by -8(0, 0) +m.\ From the model of Ref. [lOj, which 
uses a model gluon propagator that is the sum of a momentum-space delta function and 
the one-loop renormalisation group result together with a dressed quark-gluon vertex, one 
infers a value of c\ = (c u — Cd)/2 = 0.006 for m u = 5 MeV and rrid = 10 MeV [by matching 



6 



c(m) to the m dependence of the solutions at q = and using the C = 500, r = e, Nf = 4 



results from Ref. ||10|| 1. This is an indication of the size of the isovector component of DCSB, 
which is seen to be a 1% effect in this case. 

The mixing amplitude obtained with this parameter set is given by the dashed curve in 
Fig. 1: less than 25% of the mass-shell mixing strength is accounted for. Of course, the 
above arguments are simply a guide. If C\ is treated as a free parameter a fit to the mass- 
shell mixing amplitude requires c\ = 0.009, which yields the momentum dependence shown 
by the solid line in Fig. 1. Clearly, the quark loop mixing mechanism is extremely sensitive 
to the small isovector component of DCSB. 

Qualitatively, the momentum dependence of the mixing amplitude is insensitive to the 
form of the model Bethe-Salpeter amplitudes at the vertices. The insert in Fig. 1 illustrates 
this by comparing the previously mentioned result [solid line] from the simple gaussian 
amplitude to the result [dash-dot line] obtained using the dominant p-amplitude obtained in 
a recent Bethe-Salpeter calculation ||18|| . Both representations for the on- mass-shell meson 
vertex give essentially identical results. 

In Fig. 2 we isolate the part of the off-mass- shell behaviour of the mixing amplitude that 
is produced by the momentum dependence of the function Z(P 2 ). This is a measure of the 
extent to which the diagonal meson propagators depart from the standard point form due 
to their qq substructure. The solid curve is the result previously shown and includes the off- 
shell contributions from the propagators, while the dashed curve has Z(P 2 ) held fixed at the 
mass-shell value to remove them. Clearly this contribution is significant for P 2 > 1 GeV 2 . 

6. CSB potential. In Fig. 3 we display the CSB NN potential V pui {r) calculated from 
the standard momentum space form [p]] 

v (p s _ 9 P N 9uuN F pN (P 2 ) II T (P 2 ) F uN (P 2 ) 
M ' (P 2 + m 2 )(P 2 + m 2 J ' 1 ' 



In accordance with the Bonn meson exchange NN interaction model [22], monopole form- 
factors for pNN and uNN coupling are used with a common range of 1.5 GeV, and the 
coupling constants are taken to be g 2 N /4ii = 0.41 and g 2 N /4n = 10.6. The mixing ampli- 
tudes in Fig. 2 give rise to the potentials shown by the solid and dotted curves in Fig. 
3. The small difference between them indicates that the off-shell structure of the diagonal 
propagators of the composite vector mesons does not make a signficant contribution to the 
potential. This is because Z(P 2 ) departs significantly from its mass-shell value only for 
momenta that are well suppressed by the meson-nucleon form factors. Variations of ±20% 
in the form-factor range do not change this conclusion. The dashed curve is the result that 
follows from the mass shell assumption, Ut(P 2 ) ~ My). 

The typical quark loop result for the nuclear CSB potential is very much weaker than 
the potential produced by the mass-shell assumption because the mixing amplitude near 
P 2 pa is about four times smaller than the empirical mass-shell value. This is a general 
consequence of a node at very low momentum and an almost linear connection with the 
mass-shell point. 

There are uncertainties in QCD-models at this scale but our investigations have shown 
no indications that the main conclusion above is affected. For example, the off-mass-shell 
behaviour of the meson-quark vertex is not guaranteed to be represented well by a single 
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mass-shell Bethe-Salpeter amplitude. In this regard, we have employed an estimate for the 
relevant P 2 > extrapolation of the gaussian representation of T(q;P) described earlier. 
This can be defined within the GCM by the variational approach |15] to the eigenvalue 
problem based on fl5|). The dominant effect is incorporated by a momentum-dependent 
range a(P 2 ) for the q 2 behaviour that increases with increasing spacelike momentum. The 
net result is an increased spacelike suppression of the mixing amplitude with a typically 
smaller value at P 2 « 0. 

7. Quark loop generates weak mixing. The quark loop mechanism will, in general, 
produce a weak mixing at small spacelike momenta. This can be seen from an analysis of 
the loop integral in @ and (|10J) at first order in the isovector components of a v and a s for 
P 2 > 0. 

From (0), both the isoscalar and isovector <j v are always positive. Thus the vector con- 
tribution to Ut{P 2 ) is positive for sufficiently large P 2 > 0. From flTTp, the isoscalar a s is 
always positive, while the isovector component has a small positive contribution from the 
first term, and a dominant negative contribution from the second term. Thus the scalar 
contribution to Ut{P 2 ) is also positive. Since the vector contribution is dominant, there 
will be a change in sign at small momenta driven mostly by the sign change in the coeffi- 
cient of the vector part of the integrand in fllPp. The position of the node is approximately 
Pq ~< 4s(l + 2u 2 )/3 > with the expectation value estimated from the dominant behaviour 
of the remaining integrand in flTUD. Using only the typical momentum scales that charac- 
terise the propagator and the vertex, we obtain Pq ps 0.5 GeV 2 which is in accord with the 
numerical results. 

On the time-like side of the node, the linear factor of P 2 in (|10D represents the behaviour 
well. 

These observations hold equally well if a free constituent quark propagator is used. 

8. Summary and conclusions. We have calculated the quark loop contribution to the 
p°-u) mixing amplitude in a way that retains, for the first time, a number of self-consistent 
features that follow from use of a QCD-based model field theory. These include DSE model 
guidance for the quark propagator at the time-like and complex momenta that arise, calcula- 
tion of the strength and range of the vertex function within the model, and the contribution 
to the off-shell behaviour arising from the composite nature of the diagonal meson propaga- 
tors. 

The mass-shell value of the mixing amplitude is found to be very sensitive to the isovector 
component of DCSB. For the behaviour at spacelike momenta relevant to the CSB compo- 
nent of the NN interaction, none of the new elements we include alter the conclusion that 
the quark loop mechanism alone generates an insignificant CSB potential. 

This failure suggests that other isospin-symmetry breaking intermediate states, such 
as the virtual pion-loop process: p — > nn — ► to through a G-parity violating umr vertex, 
may contribute significantly to Ut{P 2 )- In this connection, and in the context of the field 
theoretical model considered herein, the unit vertex occurs as a quark loop which is non- 
zero for m u ^ rrid and a recent estimate |L6| indicates that g^nn/gprnr ~ ^(-My)/™^. It 



is therefore possible that the pion loop mechanism could make a significant contribution to 
p° — uj mixing and the associated NN CSB potential should be investigated. This is especially 
true since the imaginary part of such an amplitude relates to the same interference effect in 



S 



the In decay channel that is the basis of the accepted experimental value of p° — u mixing. 
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FIGURES 



FIG. 1. The quark loop contribution to the p° — u> mixing amplitude. Here P 2 > represents 
spacelike momentum. The solid curve is obtained with the quark propagator parameter c\ adjusted 
to fit the experimental mixing value; the dashed curve uses the value inferred from a realistic DSE 
solution; the dotted curve is produced by the common assumption that a current mass provides 
a constant shift to the scalar self-energy function. The insert demonstrates that the mass-shell 
vector Bethe-Salpeter amplitude does not introduce significant model dependence. 

FIG. 2. The p° — u> mixing amplitude with (solid curve) and without (dashed curve) the 
off-mass-shell contribution from the calculated diagonal propagators of the qq mesons. 

FIG. 3. The p° — to mixing contribution to the CSB NN potential form factor. The dashed curve 
is obtained under the assumption that the mixing amplitude is a constant given by the mass-shell 
value. The quark loop calculation is shown by the solid curve when the off-mass-shell contribution 
from the diagonal meson propagators is included, and by the dotted curve when it is not. 
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